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ABSTRACT 

Warm Neptune- and sub-Neptune-sized exoplanets in orbits smaller than Mercury’s are thought to have 
experienced extensive atmospheric evolution. Here we propose that a potential outcome of this atmospheric 
evolution is the formation of helium-dominated atmospheres. The hydrodynamic escape rates of Neptune- and 
sub-Neptune-sized exoplanets are comparable to the diffusion-limited escape rate of hydrogen, and therefore 
the escape is heavily affected by diffusive separation between hydrogen and helium. A helium atmosphere can 
thus be formed - from a primordial hydrogen-helium atmosphere - via atmospheric hydrodynamic escape from 
the planet. The helium atmosphere has very different abundances of major carbon and oxygen species from 
those of a hydrogen atmosphere, leading to distinctive transmission and thermal emission spectral features. 

In particular, the hypothesis of a helium-dominated atmosphere can explain the thermal emission spectrum of 
GJ 436 b, a warm Neptune-sized exoplanet, while also consistent with the transmission spectrum. This model 
atmosphere contains trace amounts of hydrogen, carbon, and oxygen, with the predominance of CO over CH4 
as the main form of carbon. With our atmospheric evolution model, we find that if the mass of the initial 
atmosphere envelope is 10“^ planetary mass, hydrodynamic escape can reduce the hydrogen abundance in the 
atmosphere by several orders of magnitude in ^ 10 billion years. Observations of exoplanet transits may thus 
detect signatures of helium atmospheres and probe the evolutionary history of small exoplanets. 

Subject headings: radiative transfer — atmospheric effects — planetary systems — techniques: spectroscopic 
— planets and satellites: individual (GJ 436 b) 


1. INTRODUCTION 

Recent exoplanet surveys have discovered warm Neptune- 
and sub-Neptune-si zed planets in tightly bound orbits around 
their pa rent stars (iButler et al.l l2004t lEressin et alJ 120131: 
[HowardI l2013h . Were these planets in our Solar System, 
their orbits would be interior to that of Mercury. The 
atmospheres of Neptune-sized exoplanets are assumed to 
be similar to those of giant planets in o ur Solar System, 
i.e., dominated by hydrogen and helium (iRoge rs & Se agerl 
2010t iNettelmann et al.ll2010l: iMadhusudhan & SeagedUOl 1 : 
Line et al.ll2014ft . Howeyer. this assumption may not be yalid 
for warm extrasolar Neptunes in close-in orbits. Due to stellar 
irradiation, these pla nets may haye experienced significant a t- 
mospheric eyolution (iLopez et al.ll2012HOwen & Wull201^ . 

Here we consider if some Neptune- and sub-Neptune-sized 
exoplanets can haye their atmospheres depleted in hydro¬ 
gen but abundant in helium, due to extensiye atmosphere 
loss. Eigure [T] illustrates a potential eyolution scenario for 
an exoplanet that leads to an atmosphere dominated by he¬ 
lium. The planet starts with a hydrogen-helium atmosphere 
accreted from the planet-forming nebula. This initial atmo¬ 
sphere may eyaporate and reduce its mass by accretion heat¬ 
ing, impact erosion, stellar winds, and X-ray and extreme ul- 
trayiolet (EUV) radiation from the parent star. Subsequently, 
the planet may conti nue to experience atmospheric loss due 


to stellar irradiation ( 

Lammer et aLII2003t lErkaev et al.l 

2007; 

Baraffe et al.l 120081 

Lonez et al.l 120121 lOwen & Wul 

2013; 

Inamdar & Schlichtingl 120141). If the escape rate is 

very 


high, the planet may lose both hydrogen and helium com- 
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pletely. Howeyer, if the escape rate is comparable to the 
diffusion-limited escape rate of hydrogen, the escaping flow 
may be highly enriched in hydrogen that is 4 times lighter 
than helium. We show later in this paper that Neptune- and 
sub-Neptune-sized exoplanets undergoing transonic hydrody¬ 
namic escape would generally haye an escape rate comparable 
to the diffusion-limited escape rate of hydrogen, and the frac¬ 
tion of hydrogen in their atmospheres would decrease as the 
planets eyolye. Our model suggests that if the initial mass of 
the atmosphere is 10“^ planetary mass, the remaining atmo¬ 
sphere would become helium-dominated within ~ 10 billion 
years. 

This concept of fossil helium atmosphere can be important 
because sub-Neptune-sized planets are ubi quitous in our in¬ 
terstellar neighborhood (e.g. lHowardll2013l) , and because he¬ 
lium atmospheres would haye yery different molecular com¬ 
positions from hydrogen atmospheres, leading to distinctiye 
spectral features in transmission and thermal emission. The 
chemical compositions and spectral features of the atmo¬ 
spheres on super Earths and mini Neptunes are highly sen- 
sitiye to the elemental abundances of the atmospheres, in par¬ 
ticular to the hydrogen abundance and the carbon to oxygen 
(C/O') ratio (iMoses et ahllMT^ iHu & Seageill2014ft . Com¬ 
pared to other non-H2-dominated exoplanet atmospheres, a 
helium atmosphere has a more extended scale height, and 
presents a better opportunity to characterize highly eyolyed 
exoplanet atmospheres yia transmission spectroscopy. 

The paper is organized as follows. We describe our models 
in § 2. We deyelop an atmosphere eyolution model to study 
the conditions to form helium atmospheres on Neptune- and 
sub-Neptune-sized exoplanets, and upgrade an existing atmo- 
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spheric chemistry and radiative transfer m odel for exoplanets 
(iHu et al.ll2012ll^l3l:lHu & Seageill2014l) to include the ca¬ 
pability of treating helium atmospheres. We present a gen¬ 
eral result of the fractionation between hydrogen and helium 
for exoplanets undergoing transonic hydrodynamic escape in 
§ 3. We then apply our hypothesis to the Neptune-sized exo¬ 
planet GJ 436 b in § 4, and propose that the planet can have 
a helium atmosphere, which explains its puzzling emission 
features observed with Spitzer. We discuss how transmission 
spectroscopy of exoplanets can distinguish highly evolved, 
helium-dominated atmospheres from hydrogen atmospheres 
in § 5, and we conclude in § 6 . 


2. MODEL 

2.1. Atmosphere Evolution Model 

We have created an atmosphere evolution model to inves¬ 
tigate whether a primordial hydrogen-helium atmosphere on 
short-period exoplanets can evolve to a helium-dominated 
one. We trace the atmosphere mass and elemental abundances 
from the present time to the past, focusing on the effects of at¬ 
mosphere escape. This way, our model can be used to test if 
a helium atmosphere is a plausible evolutionary outcome for 
an exoplanet. We truncate the model at an age of 0.1 bil¬ 
lion years, because the processes dominating the period less 
than 0.1 billion years after the planet formed (e.g., accretion 
heating, impact erosion, rapid photoevaporation) do not frac¬ 
tionate hydrogen versus helium. To simplify the calculation, 
we assume that the core to be unchanged in mass, size, and 
composition in the evolution history, and that the atmosphere 
is isothermal. The cor e could cool down a nd shrink by up to 
1 R 0 over the history (iBaraffe et alj|200^ . We have verified 
that the evolution history of the atmospheric hydrogen abun¬ 
dance is not qualitatively sensitive to these assumptions. 

A key element of the evolution model is escape of 
hydrogen and helium from the envelope. For irra¬ 
diated exoplanets, the escap e is mainly driven b y the 
stella r EUV radiation (e.g. iLammer et al.l 120031: lYelle 


2004t [Lecave lier des Etan gs et akl 20041; Tian et alj 


Garcia Munozn2007r Murrav-Clav et al.( l2009t iGuol _ _ .. 
Koskinen et al.ll20T^ Lammer et al.1120 1 jfe . and the escape 


rate can be estimated by an energy-limited escape formula that 
relates the EUV irradiation to adiabatic expansion and cool¬ 
ing of the outflow, with an efficiency ( 77 ) that umbrellas the 
unde rlying complex ene rgetic and hydrodynamic processes 
(e.g. lErkaev et al.ll2007h . Typical values for rj vary from a 
few percent to a few tens percent and may change during 
the evolution dependi ng on the composition and the ioniza¬ 
tion d egree of the flow (iMurrav-Clav et al.l2009HOwen & 
l2Q13h . It has also been suggested that X-ray irradiation may 
drive escape very early in the history when the ste l lar X- 
ray luminosity is high (ICecchi-Pestellini et al.1120061120091 : 
lOwen & JacksonI l2012h . However, X-ray-driven escape is 
less important for our purpose than EUV-driven escape be¬ 
cause much of the X-ray en ergy deposition is reradi ated away 
and does not drive escape (lOwen & .Tacksonll2012h . and be¬ 
cause the X-ray flux is one-order-of magnitude smaller than 
the EUV flux after 1 billion years since formation (see be¬ 
low). 

The energy-limited escape rate is 


_ lEuvna^rl 

~ AKd^GM„ ’ 


( 1 ) 


in which Leuv is the luminosity of the star in EUV, t] is the 


heating efficiency, or the fraction of absorbed EUV energy 
that drives escape, a is the ratio between the planet’s EUV 
absorbing radius and the planetary radius, is the planetary 
radius, Mp is the planetary mass, d is the semi-major axis, 
G is the gravitational constant, and K is the pot ential energy 
reduc tion factor due to the Roche lobe effect (lErkaev et all 
[2003). Eor convenience we define Q as the incident EUV 
heating rate on a planet, i.e. 


Le^ 
^ Ad^ ’ 


( 2 ) 


and 2net = pQ as the net EUV heating rate that drives escape. 

Eor Leuv, we can adopt the following empirical estimates 
for stellar EUV fluxes based on X-ray observations and coro¬ 
nal modeling of 80 stars with s pectral types spanning from M 
to E (ISanz-Forcada et aT]l201 Ih 


log(LEuv) = 22.12-1.241og(r), (3) 

log(Lx-ray) = 21.28 - 1.44 log(r), (4) 


where the luminosity has a unit of J s“' and the age (r) has 
a unit of billion year. This estimate of the EUV luminosity 
has an uncertainty of an order of magnitude for any individ¬ 
ual star; however, we will see later that our result does not 
sensitively rely on the EUV luminosity. Equations (|2l and (|4|i 
indicate that the X-ray luminosity can only be a minor contri¬ 
bution to the escape energy for the bulk part of the evolution. 

The planetary radius is defined to be the radius of the ho¬ 
mopause, i.e., the level above which the binary diffusion co¬ 
efficient between hydrogen and helium is greater than the 
eddy diffusion coefficient. Such a definition is only nomi¬ 
nal because significant hydrogen photochemistry may occur, 
and the mixin g ratio of H is not a constan t below the ho¬ 
mopause (e.g. Li ang et al.ll2003l lYellell2()()4t iHu et al.l 1201 21 : 
iKoskinen et al.ll2013l) . For a n eddy diffusion coeffic ient rang¬ 

ing in 10® ~ lO" cm^ s“' (iParmentier et al.l[20131) . the ho¬ 
mopause is at 10“^ ^ 10“^ Pascal for a temperature of 600 K, 
and 10“' ^ 10“^ Pascal for a temperature of 10"' K. 

Most of the EUV irradiation is deposited at the r = 1 pres¬ 
sure level (lYung & Demoreiri999l) . 


1 

Pi = nikT ^- kT 

aH 


mug 

a 


(5) 


where H is the scale height, ma is the mass of hydrogen, g is 
the gravitational acceleration, and cr = 2 ~ 6 x 10 “'^ cm^ is the 
cross section of hydrogen for the EUV irradiation in 13.6-20 
eV (ISpitzeiill978l) . For GJ 436 b (g ~ 10 m s“^), the pressure 
level is 3 ^ 8 X 10“^ Pascal. Therefore, the EUV absorbing 
altitude is up to 9 scale heights above the homopause, which 
corresponds to less than 0.1 planetary radii. In other words, 
a < 1.1. We therefore assume a = 1 in this work, effectively 
assuming the uncertainty in in Equation ([T]i is absorbed in 
a wide range of rj. The parameter a being clos e to unity is 
also verified by full thermosphere models (e. g. lYellel 120041 : 
iMurrav-Clav et al1l2009t iKoskinen et al.ll2013l) . We note that 
this assumption does not imply EUV heating occurs at the 
homopause, and the parameter a can be greater for smaller 
planets (e.g. ILammer et al1l2013l) . 

In order to apply Equation ([1]) to calculate the escape rate, 
we need to determine the regime for the hydrodynamic out¬ 
flow. Hydrodynamic outflow from planetary atmospheres 
powered by external heating has recently been calculated by 
Monte Carlo direct simulations, and the advantage of these 
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Initial state 



After rapid evaporation 
and impact removal 



Asymptotic state 



t ~ several Gyr 


Figure 1. An evolution scenario for an in'adiated Neptune- or sub-Neptune-sized exoplanet that produces a helium-dominated atmosphere from a primordial 
hydrogen-helium atmosphere. Rapid evaporation could last for ~ 0.1 Gyr, depending on the thermal history of the planet and the evolution of the stellar X-ray 
and EUV radiation, and the subsequent thermal escape modulated by diffusive separation bet\veen hydrogen and helium could last for multiple billion years to 
form a helium-dominated atmosphere. The required time depends on the mass of the initial atmosphere. 


simulations is that they no longer require outer boundary con¬ 
ditions and that they self-consistently determine the location 
of the Knudsen layer be tween the static and the hydrodynamic 
parts of the atmosph e re (iTu cker et a l.l2012l:lErwin et al J2013l : 
iJohnson et alJ l2013t IVolkov & JohnsonI 120131) . These sim- 
ulations find the energy-limited escape rate to be a good 
approximation of the escape rate when the escape flow is 
subsonic, but not when the escape flow becomes transonic. 
The minimum heati ng rate to drive a transonic outflow is 
jJohnson et alj|20l'3l) 


Qnet > Qc 47rr* 


7 


C r O'rfCllft 


2 t/(r,) 


U(rn), 


(6) 


where is the critical heating rate, r* is the sonic point ra¬ 
dius, 7 is the heat capacity ratio, Cc<Jc is the collisional cross 
section, Kn„t is the maximum Knudsen number for the flow to 
be in the continuum regime, U is the gravitational energy of 
the escaping particle, and m is the mass of the particle. For an 
atomic flow, we adopt 7 = 5/3, CcCTc = 5 x 10“^° m^, Kn„, ^ 1 
(valid when the heat is primarily absorbed over a broad range 
of radius below r*), and r, ^ tq. Note that r* can be greater 
than rp, but Qc only weakly depends on its specific value as 

In the transonic escape regime, the escape rate no longer 
increases with the energy put; instead, incremental energy 
input would be taken awa y as thermal and translational en¬ 
ergy (IJohnson et al.ll20f3h . Therefore, Equation ([T]i can sig¬ 
nificantly overestimate the escape rate for transonic flow. We 
instead estimate the escape rate as 


$=/.$EL, (7) 

where fr, adopted from IJohnson et al.l (120131) . is the escape 
rate reduction factor due to the flow being transonic, fr is 
a function of QaeilQc' when QnmlQc < f /r ~ 1 and we 
return to the conventional energy-limited escape; however, 
when QnmlQc > 1 . /r < 1 and it decreases nearly as inverse- 
proportionally to a greater QaeilQc- 
Equation (|7]i allows us to trace the history of atmosphere 
loss in terms of the total mass loss flux. Next, in order to trace 
the effect of atmosphere loss on the atmosphere’s elemental 
abundance, we need to partition the flux into the mass flux of 


escaping hydrogen and the mass flux of escaping helium, i.e., 

$ = =47rr^((/H»JH+'/He»JHe), (8) 

where (/h and (/hb denote the number fluxes, and mn and 
mHe denote the masses of hydrogen and helium atom, respec¬ 
tively. The escaping flow should be made of atomic hydrogen 
and helium, because photochemistry models show that the 
dominant form of hydrogen in the upper atmosphere should 
be atomic hydrogen, due to photod issociation of hydrogen 
molecules catalyzed by water vapor (iLiang et al.ll2003l : lYelld 

IliollMoses et akll^llHu et al.ll2Q12l). 

The physical reason for different fl>H and <i)He is that the he¬ 
lium atom is 4 times more massive than hydrogen and that 
both fluids are subject to the gravity of the planet. When $ 
is large, the coupling between the two components is strong 
and therefore the two fluxes are approximately proportional 
to the mixing ratios of the gases at the homopause, which re¬ 
sults in little fractionation. When $ is small enough helium 
will not escape and hydrogen escape will be in the so-called 
“diffusion-limited” escape regime. The nature of mass frac¬ 
tionation in hydrodynamic escape has bee n investigated thor¬ 
oughly with hydrodynam ic formulati ons (IZahnle & Kastingl 
[T^lHunten'S^fT^IZahnle et allfTom 

Two complications arise as we apply the classical calcula¬ 
tion of mass fractionation to irradiated exoplanets. We have 
shown earlier that the homopause can be quite close to the 
EUV absorption radius, which means that (1) the tempera¬ 
ture of the homopause can be much higher than the planetary 
equilibrium temperature due to EUV heating; and (2) the out¬ 
flow can be partially ionized and the ion-neutral interactions 
should be considered in assessing whether escaping hydrogen 
can drag helium. 

The rate of momentu m exchange between hy drogen and he¬ 
lium can be written as (ISchunk & Nagvlll98d) 


(jMne 

5t 


kT ^ ^ mn+v 

«He(MH - MHejWH 7— + «He(MH+ “ MHe)«H+-, 

b WHe 


(9) 


where n denotes the number density, u denotes the veloc¬ 
ity, b = 1.04 X cm“' s“* is th e binary diffusion 

coeffi cient between hydrogen and helium (iMason & Marrero! 
[I9^ . and V is the ion-neutral m omentum transfer collision 
frequency (ISchunk & NagylflOSOl) . The first term is the mo¬ 
mentum exchange between neutral hydrogen and helium, and 




















































4 


Hu et al. 



Temperature [K] 


Figure 2. Comparison between the H-He coupling and the H'''-He coupling. 
The H-He coupling ffle term kT jb in Eg. [3 is calculated from the binary 
diffusion coefficient iN^son & Marrer3ll97^ . and the H'''-He coupling (the 
term mn+v/nHe in Eq.[9) is calculated from the ion-neutral momentum trans¬ 
fer collision frequency <Schunk & Nag’vll98(]n . 


(IZahnle et al.lll99dl) . This equation demonstrates that when 
(/)H 3> Xii(j)uL, fractionation between H and He would be min¬ 
imal, but when 0 h ^ </'He —S’ 0, we return to the classi¬ 

cal formula for diffusion-limited escape. Near this limit, both 
H and He escape, so the flow is not exactly diffusion-limited 
escape, but significant fractionation between H and He can 
still occur. With Equations (0 and fOT l we determine the es¬ 
cape flux of hydrogen and helium. The solution is 

If $ < (/)DL^H'«H47rr^, 

4>h = 4>, 4>He = 0; (14) 

If $ > (/)DLi^H'«H47rr^, 

^mnXy, + 4>DhmHmiieXiiXii^4Trrl 

4*h =- 

niYlX^ + rWHe-^He 

4'He= - -■ (15) 

For each time step, the mass of the atmosphere is updated 
by adding the escape fluxes of hydrogen and helium, and 
the composition and the mean molecular mass of the atmo¬ 
sphere are updated accordingly. Then, the new mass and mean 
molecular mass are used to determine the new thickness of the 
atmosphere and the new radius of the planet. 


the second term is the momentum exchange between ionized 
hydrogen and helium. Helium may also be partially ionized, 
but t he ionization fraction of helium is low up to large radius 
(e.g. iKoskinen et ani2013h . The relative efficiency between 
the two exchange processes is shown in Figure |2] We see 
that ionized hydrogen or heated hydrogen can better drag he¬ 
lium than cold neutral hydrogen, both by up to a factor of 2. 
Assuming mh+ = mh, the rate of momentum exchange can be 
written as 


JMhe 

5t 


«He(MH - MHe)(«H + «H+) 


kT 

V' 


( 10 ) 


where b' is the effective binary diffusion coefficient, defined 
as 


kT 

V 


(1 -x)—+x 
b 


oth+v 

MHe 


( 11 ) 


where x is the ionization fraction. As such, the diffusion- 
limited escape rate is 


0DL = 


rjkT 


( 12 ) 


where T is the temperature of the homopause. 

The fact that the efficiencies of the neutral interaction and 
the ion-neutral interaction are similar and only weakly depend 
on temperature (Figure |2]i allows us to choose a representa¬ 
tive temperature and ionization fraction. As a conservative 
estimate for the mass fractionation effect, we use a high tem¬ 
perature of 10"^ K, and an ionization fraction of 0.1, which 
yields b' = 8.0 x 10^° cm“* s“* based on Equation (fTTl l 

The diffusion-limited escape rate determines the difference 
between the two escaping components, as 


0He 

^He 



(13) 


This simple relation is derived assuming subsonic flow but is 
proved to be a close approximation to a transonic flow as well 


2 .2. Atmosphere Chemistry Model 

Atmospheric escape and evolution can have large effects 
on the atmospheric chemical composition and then a planet’s 
transmission and thermal emission spectrum, by depleting 
hydrogen over a long period. The fact that vertical mixing 
(’^years) is much faster than atmospheric evolution (^billion 
years) allows us to treat the atmospheric chemistry calcula¬ 
tions as “snapshots in the evolutionary history”, with the in¬ 
gredients of the atmosphere (i.e., the elemental abundances) 
controlled by the evolution. 

We use a one-dimensional photochemistry- 
thermochemistry kinetic- transport model (IHu et al.l 120121 
l2013t [Hu & Seag^l2014ll to explore the chemical state and 
spectral characteristics of a potential helium atmosphere. 
The photochemistry-thermochemistry kinetic-transport 
atmosphere model computes the atmospheric composition 
and the opacities self-consistently for an exoplanet, with the 
atmospheric elemental abundance, the efficiency of vertical 
mixing, and the internal heat flux as the input parameters. 
The model also computes the temperature profiles from 
the opacities using the grey-atmosphere approximation. 
Our atmosphere chemistry model contains two steps of 
computation: the first step is to compute the atmospheric 
composition at thermochemical equilibrium self-consistently 
with the temperature profile; and the second step is to use the 
result of the first step as the initial condition to simulate the 
effects of vertical mixing and photochemical processes on 
the atmospheric composition, with the temperature-pressure 
profile adjusted accordingly. After the photochemistry- 
thermochemistry simulation converges to a steady state, we 
compute synthesized spectra of the modeled exoplanet’s 
atmospheric transmission and thermal emission with a 
line-by-line method. A unique feature in our photochemistry- 
thermochemistry model is that the model does not require 
specification of the main component of the atmosphere (nor 
the mean molecular mass) and instead the model takes the 
elemental abundances as the input parameters, making our 
model ideal for exploring the composition of hydrogen- 
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versus helium-dominated atmospheres. 

Our general atmosphere chemistry mo del and i t s val¬ 
idation have been descr ibed in detail (iHu et al.l 120121 
120131 iHu & Seageri 120141) . To properly treat a helium- 
dominated atmosphere, we have made the following up¬ 
grade to the general model. We have included th e H 2 - 
He collision-induced absorption (iBorvsow et al.l[T988h in our 
calculations of the opacities, and used the heat capacity 
of helium recommend by the NIST Chemistry Webbook 
(http://webbook.nist.gov/chemistry/) for calculation of the 
adiabatic lapse rate of helium-dominated atmospheres. We 
have treated helium as a chemically inert gas in our model, 
and included helium in calculation of the mean molecular 
mass. 



10 10 10 10 10 10 10 
Q^Q, i.e. the mininum heating efficiency to drive transonic escape 


3. FRACTIONATION BETWEEN HYDROGEN AND HELIUM 

We Study whether hydrodynamic escape of hydrogen can 
fractionate hydrogen and helium in an exoplanet atmosphere. 
First, most of the detected exoplanets receive enough EUV 
radiation to have transonic hydrodynamic escape. Figure |3] 
shows that a majority of the detected exoplanets with their 
masses and radii measured would undergo transonic hydrody¬ 
namic escape for a moderate heating efficiency of ~ 10% at an 
age of 10 billion years. When they are younger and the stellar 
EUV luminosities are higher, the planets are even more likely 
to experience transonic escape. This is reasonable because 
most of the planets have short orbital periods and locate close 
to their parent stars. Therefore, the transonic regime must be 
considered when estimating the hydrodynamic escape rate of 
the short-period exoplanets. 

Second, hydrodynamic escape from many Neptune- and 
sub-Neptune-sized exoplanets significantly fractionates hy¬ 
drogen versus helium in their atmospheres. Assuming tran¬ 
sonic hydrodynamic escape, we can replace Qnet by in 
the energy-limited escape rate calculation for an estimate of 
the upper-limited of the escape mass flux, because increas¬ 
ing the stellar EUV flux would n o longer increase the escape 
mass flux for transonic outflow dJohnson et al.ll20f3h . Note 
that this upper limit is independent from the stellar irradia¬ 
tion, but depends on the mass and the radius of the planet, as 
well as the composition of the atmosphere. We then partition 
this mass flux to the escape fluxes of hydrogen and helium, 
using Equations (fl4l and ffSl l. and calculate the fractionation 
factor between hydrogen and helium, x^, which is defined as 


Y2 = 


0He/0H 

Vhc/Vh 


(16) 


A 2 depends on the instantaneous volume mixing ratio between 
hydrogen and helium, and takes a value between 0 and 1. X 2 = 
1 means no fractionation, and X 2 = 0 means the highest degree 
of fractionation (i.e., no helium escaping with hydrogen). We 
show the results in Figured 

Most of the Neptune- and sub-Neptune-sized planets have 
a small fractionation factor between hydrogen and helium. 
This means that hydrodynamic escape from their atmospheres 
would deplete hydrogen and enrich helium. Since the as¬ 
sumed escape rate is an upper limit, we are conservative in 
estimating this fractionation effect. For a smaller escape rate, 
the fractionation effect would be even more significant. We 
also note that the fractionation effect becomes milder as the 
hydrogen mixing ratio decreases. However, the asymptotic 
value of X 2 is less than unity (Figure|4|). This implies that with 
unlimited time and no hydrogen source, the process can in 


Figure 3. The minimum heating efficiency to drive transonic hydrodynamic 
escape from detected exoplanets, shown as a histogram. The exoplanets hav¬ 
ing their masses and radii measured more precisely than 3 - (T are counted, 
and the total number of planets are 254 (http://exoplanet.eu I. We use Equa¬ 
tion (D to estimate the incident EUV heating rate, and assume a conservative 
10-billion-year-old stellar EUV luminosity from Equation (3). We use Equa¬ 
tion (B to estimate the critical heating rate for transonic escape, assuming 
that the sonic point radius is close to the planetary radius, and that the mean 
molecular mass is 2 (i.e. 67% H and 33% He). 


principle remove hydrogen in the atmosphere to an arbitrarily 
low abundance. 

Therefore, fractionation between hydrogen and helium is 
significant for hydrodynamic escape from many exoplanets. 
This does not imply these planets have helium atmospheres, 
however, as additional conditions have to be met. Our at¬ 
mosphere evolution model serves to explore the conditions 
to form helium atmospheres for specific planets. Neptune- 
and sub-Neptune-sized planets can have small hydrogen and 
helium envelopes, and the fractionation factors for them are 
substantially smaller than unity. Note that the fractionation 
factors for them are not zero, so they are technically not in 
the regime of diffusion-limited escape. Rather, hydrodynamic 
escape from these planets is highly affected by the diffusion 
separation between hydrogen and helium. We suggest that 
fractionation between hydrogen and helium is an important 
process for hydrogen-helium atmospheres of short-period ex¬ 
oplanets undergoing hydrodynamic escape. The effect of this 
fractionation, as later discussed for the cases of GJ 436 b, de¬ 
pends on the initial mass of the atmosphere, and is likely to 
be significant for Neptune-sized and smaller exoplanets. 

4. A HELIUM ATMOSPHERE MODEL FOR GJ 436 B 

The observational characterization of atmospheric com¬ 
positions of Neptun e- and sub-Neptune-size d exoplan¬ 
ets is ongoing (e.g. IS eager & Demind 120101) . and pro¬ 
vides the opportunity to search for highly evolved exo¬ 
planet atmospheres and test our hypothesis of the forma¬ 
tion of helium atmospheres. One of the most observed 


Neptune-sized exoplanets is GJ 436 b (iGillon et al.l 

2007; 

Deming et aklUGOTl; iDemorv et al.ll2007l; lAlonso et al.l 

2008; 

Pont et alJ 20091: Caceres et alJ 20091 IBallard et al.l 

201C; 

Stevenson et all 120101: Beaulieu et al.l 1201 It Knutson et al. 

201 11 12014at ILanotte et alll2014l IMorello et al.ll2015l). 

This 


planet orbits a nearby M dwarf star and has a mass of 23.2 Mq 
and a radius of 4.22 (iTorres et al.ll2008l) . The bulk density 
of the planet implies that the planet should have a thick atmo¬ 
sphere or envelope made of hydrogen and/or helium, which 
itself has a mass of at le ast 10“^ ~ 10“^ of the planet’s mass 
(iNettelmann et al.ll2010l) . 
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Figure 4. Fractionation factor between hydrogen and helium for Neptune- 
and sub-Neptune-sized exoplanets undergoing transonic hydrodynamic es¬ 
cape. The color contour shows the fractionation factor (X 2 ), for which a 
smaller value means a more significant fractionation between hydrogen and 
helium. Neptune- and sub-Neptune-sized exoplanets having their masses and 
radii measured more precisely than 3 — cr are shown. 


4.1. Puzzling Emission Features of GJ 436 b 

The composition of GJ 436 b’s atmosphere has been 
a puzzle. A broad-band emission spectru m of the planet 
in 3G0 ftm was obtained from Spitzer (IStevenson et al.l 
l2010t iLanotte et al.l 120141) . which indicated that the 
planet’s atmosphere is poor in CH^ and ri ch in CO 
(iMadhusudhan & Sealed 1201 ll iLine et aP 120141) . The key 
spectral feature is a detected emission flux in the 3.6 pm 
band (CH 4 absorption band) and non-detection in the 4.5 pm 
band (CO absorption band). We note tha t the reported 3.6 
pm emission flu x of IL anotte et ^ (1201 4|) was signihcantly 
lower than that of iStevenson et al.l (120101) . though both values 
were derived from the same observation (see Appendix 0 
for discussion). In addition, a recent transmission spectrum 
of the planet in 1-2 pm was obtain ed from HST, and th e 
spectrum lacks H 2 O or CH 4 features (iKnutson et al.ll^l4ah . 
These observations and their interpretation lead to a theo¬ 
retical challenge: even when accounting for disequilibrium 


processes (iLine et al.ll201 ll) . atmospheric chemistry models 
uniformly predict most carbon should be in the form of 
CH 4 -not CO-in a solar-abundance hydrogen and helium 
background atmosphere. 

A solar-abundance hydrogen and helium atmosphere on 
GJ 436 b cannot produ ce a spectrum consistent with either 
datase t of the spectra (IStevenson et al.l 120101 : ILanotte et alJ 
l2014ft . because the chemical equilibrium models would pre¬ 
dict most of carbon to be in the form of CH 4 , but the 
observations indicate otherwise. Including disequilibrium 
processes like strong vertical transport and efficient pho¬ 
tochemical reactions, or exploration of temperature vari¬ 
ation in reasonabl e ranges, cannot address the “missing- 
methane” problem (IMadhusudhan & SeagedbOl UlEIne et al.l 


I 2 OI It lAgilndez et akl 20141). To fit the emission spec- 
trum of IStevenson et al.l (120101) . the mixing ratio of CH 4 
between 1 and 0.01 bar would need to below ^ 1 ppm 
based on retrieval studies allowing the compositions and 
the temperature profiles to va r y freely and indep endently 
(IMadhusudhan & Seagerl 1201 lb iMoses et al.1 120131) . The 
retrieval results of htting to the more recent dataset of 
ILanotte et al.l (l2014l) are not available, but one could imagine 
that more CH 4 could be allowed in the atmosphere because of 
th e smaller 3.6-rtm ec lipse depth. 

iMoses et akl (1201 3l) presented a series of models for the at¬ 
mospheric compositions and emission s pectra of GJ 436 b , 
using the PHOENIX atmospheric model dBarman et alJ2005h 
for the temperature-pressure prohles. To obtain a mixing ratio 
of methane less than 1 ppm, a metallicity 10,000 times greater 
than the solar atmosphere would be required (IMoses et alJ 
I 2 OI 3 L Figure 12). Such a 10,000x solar atmosphere would 


be CO,-, CO-, N, 


and H 20 -dominated, 


and would have a 


very high mean molecular mass. 

However, this hypothetical world cannot be GJ 436 b, 
for such an extravagant metallicity violates the constraints 
on the bulk density of the planet in question. The ra¬ 
dius of an equally massive planet with a CO-, CO 2 - and 
H 20 -dominated atmosphere would be much smaller than 
the measured radius of GJ 436 b. Detailed models of the 
interior structures of Neptune-sized exo planets have also 


ruled out such a dense ou ter envelope ( Rogers.^^eagei 
120 lot Nettelmann et al.l201^ . Inparticular. lNettelmann et^ 


(I 2 OIC ) calculated a high-temperature endmember scenario, in 


which the temperature profile was assumed to be adiabatic 
and the temperature at 1 bar was assumed to be 1300 K (se¬ 
ries G). This scenario had a temperature as high as 5500 K 
at 1000 bar, but the minimum mass of hydrogen and helium 
for this model was still at least 0 . 01 % of the planetary mass. 
A 10,000x solar outer envelope would not be permitted even 
when a fairly high internal temperature is considered. 

In addition to raising the metallicities, IMoses et al.1 (120131) 
then assumed an ad hoc temperature prohle, and proposed 
that an H 2 - and He-dominated atmosphere 300 times more 
metal-rich than the solar atmosphere could be consist e nt with 
the emission spectrum reported by IStevenson et al.1 ([2013). 
To produce a sufficiently low mixing ratio of CH 4 , the tem¬ 
perature of such an atmosphere at 0.1 bar would have to be 
~ 1300 K, greater than radiative-convective model predic- 
tions by at least 200 K. In other words, the 300x solar model of 
IMoses et al.1 (120131) likely assumed an unphysically large tem¬ 
perature at 0.1 bar. The tidal heat due to the planet’s orbital 
eccentricity (e=0.16) is deposited at a much deeper pressure 
and c annot raise the temperature at 0.1 bar (lAgundez et al.1 
I 2 OI 4 I) . When using reasonable temperature prohles, raising 
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the metallicities to the extent still allowed by the planet’s bulk 
density cannot reduce the mixing ratio of CH 4 in the observ¬ 
able part of the atmosphere to less than 1 ppm. 

We also note that some of the self-c onsistent models fo r 
300x and lOOOx solar abundances in iMoses et al.l (1201 3b . 
without assuming an ad hoc temperature profile, appear to be 
consi stent with the more recent dataset of the emission spec¬ 
trum (iLanotte et al.ll201^ . 

Alternatively, one could explain GJ 436 b’s emission spec¬ 
trum by the aforementioned scenario of atmospheric evolution 
culminating in a helium-dominated atmosphere. Once such a 
state is achieved, the scarcity of atmospheric hydrogen will 
cause the main molecular carrier of carbon to be CO rather 
than CH 4 , explaining the emission spectrum. 

4.2. Chemistry and Spectrum 

Using the atmosphere chemistry model, we explore a wide 
range of parameters that define the atmospheric composition 
of GJ 436 b, including the hydrogen elemental abundance 
of the atmosphere (by number, denoted as Jfn), the metal- 
licity of the atmosphere (denoted as Xm), the carbon to oxy¬ 
gen ratio of the atmosphere (denoted as Xq/Xq), the intrin¬ 
sic temperature (to specify the internal heat flux, denoted as 
Tint), and the eddy diffusivity (to specify the efficiency of 
vertical mixing, denoted as TGz)- In this work we consider 
only H, He, C, and O species, so that Xhc = 1 -Xh-Xm and 
JfM =Xc + Xo. Nitrogen, sulfur, and other metal species may 
also contribute to the spectrum. Our selection suffices because 
the key observational feature to explain is lack of methane ab¬ 
sorption. For each atmosphere scenario that we consider, we 
use the thermochemistry-photochemistry model to compute 
the steady-state molecular composition from 10 ^ bar to 10 “* 
bar. We verify that the lower boundary at 10^ bar is sufficient 
to maintain thermochemical equilibrium at the lowest atmo¬ 
sphere layer for each simulation. For the stellar input spec¬ 
trum, we use th e latest HST measu rement of the UV spec¬ 
trum of GJ 436 (iFrance et al.ll2013h and the NextGen simu¬ 
lated visible-wavelength spectrum (I Allard et alJll997n of an 
M star having parameters closest to those of GJ 436 (i.e., ef¬ 
fective temperature of 3600 K, surface gravity log(g)=5.0, and 
solar metallicity). 

We find that a helium-dominated atmosphere with a low 
hydrogen abundance and close-to-solar carbon and oxygen 
abundances can adeq uately fit both d a tasets of GJ 436 b’s 
emiss ion spectrum (IStevenson et al.l 120101 : iFanotte et all 
l2014b . At the same time a v ery different model also generally 
fits the more recent dataset (IFanotte et al.ll201^ . a hydrogen- 
rich atmosphere with highly super-solar carbon and oxygen 
abundances (Figure |5] and Figure |6ll. There is a strong resid¬ 
ual discrepancy between the data and both models at 16 pm, 
and it can be further reduced if the temperature profile is al¬ 
lowed to adjust freely ('AppendixlBli. 

Both the helium atmosphere (i.e., the H 2 -poor model) and 
the H 2 -rich metal-rich model ensure the predominance of CO 
oyer CH 4 as the main form of carbon, the key indication of 
the planet’s emission spectrum. Our chemical models indi¬ 
cate that the atmosphere would be depleted in methane when 
Xh < (Figure|7]l. Under this condition, methane 

would haye to compete with water for the limited resource of 
hydrogen. This condition generalizes the resul t of a preyi- 
ous a tmospheric chemistry model of this planet (iMoses et alJ 
l2013h : instead of increasing the metallicity Xm, we suggest 
decreasing Xu could also make this condition satisfied. In 



Figure 5. Comparison of model spectra of the dayside emission of GJ 436 b 
with o bseiwations. The report ed values fo r the eclipse depths are shown in 
black iStevenson et aUDOlOl) and orange iLanotte et alJlzQl^ . The filled 
circles are the model-predicted occultation depths with the Spitzer band¬ 
pass incorporated. The wavelength ranges of major molecular absolution 
bands are shown below the horizontal axis. The spectra are based on at¬ 
mosphere models with parameters tabulated in Tabled and the correspond¬ 
ing temperature-pressure profiles are shown in the inserted panel. The thin 
dashed lines show the eclipse depths if the planet emits as a black body with 
temperatures of 500-1100 K. The red line shows the best-fit model to the 
data of IStevenson et al] i201(D . which has Xh = 10^ and an internal heat 
flux. The magenta line shows a model with the same composition but with- 
out the internal heat flux. The blue line shows the best-fit model to the data of 
ILanotte et alj 420141) . and the green line shows the a model of H.,-rich metal- 
rich atmospheres (Xh > 0.3, Xm >0.1). The thermal emission measurement 
at the 3.6 p m band constrains on th e potential planetary scenarios: the value 
reported bv IStevenson et all 420101) suggests a planet h aving an HT-poor at- 
mosphere and internal heating; but the value reported bv ILanotte et alJ 420141) 
allows both the scenarios of an H 2 -poor atmosphere and the scenarios of an 
H 2 -rich metal-rich atmosphere. 


addition, the weak absorption features at the 5.8-//m and 8.0- 
bands seen in both datasets suggest that the atmosphere 
should still have H 2 O and CH 4 , and lead to a lower bound of 
Xu (Figure I 2 I 1 . 

Both the He atmosphere and the H 2 -rich atmosphere sce¬ 
narios require an aerosol layer at 1 mbar pressure to gener¬ 
ate the obseryed featureless transmission spectrum (Figure[^. 
This is because a hydrogen and/or helium atmosphere, with a 
mean molecular mass of 2 - 4, would produce large absorp¬ 
tion features in transmission if it does not haye aerosols. We 
find that a thin aerosol layer haying a yertical optical depth 
less than 0.1 at mid-infrared wayelengths could produce the 
featureless transmission spectrum obseryed at 1-2 pm, and 
also consistent with the broadband transit depths at 3.6, 4.5, 
and 8.0 pm (Figure [^i. Such an aerosol layer would haye 
little impact to the emission spectra. A thicker aerosol layer 
is possible and consistent with the transmission spectrum. We 
haye assumed that the atmosphere to be free of thick clouds in 
the calculation of thermal emission, because the atmosphere 
would be too h ot to form water clou ds and too cold to form 
silicate clouds (ISudarsky et al.ll^03h . Eyen if thick clouds 
were present, the detected emission flux in the 3.6 pm band 
and non-detection in the 4.5 pm band still require the main 
molecular carrier of carbon to be CO and not CH 4 . 

The large emission flux at 3.6 ^m reported by 
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Table 1 

The parameters and the modeled compositions for the atmospheric models shown in Figure[5] 


Parameter 

Molecular Composition at 0.01-1 bar 

xV« 

Color 

Xh 


Xc/Xo 

Tint (K) 

He 

H 2 

H,0 

CH 4 

CO 

CO, 



I X 10-^ 

1 X 10-2 

0.95 

60 

0.9951 

4.8 X 10-5 

2.5 X 10-5 

3.1 X 10-'" 

4.6 X 10-5 

2.5 X 10^ 

2.5^ 

r 

3 X 10-5 

1 X 10-5 

0.90 

20 

0.9995 

1.4 X 10-5 

1.3 X 10-5 

8.8 X lO-** 

4.2 X 10^ 

5.1 X 10-5 

7.1“ 

m 

1 X 10-5 

1 X 10-5 

0.9997 

60 

0.9990 

4.8 X 10-* 

2.8 X 10-5 

6.5 X 10-5 

4.9 X 10^ 

3.9 X 10-5 

2.1^ 

b 

3 X I0-‘ 

1 X 10-' 

0.70 

60 

0.76 

0.17 

1.8 X 10-2 

4.1 X 10^ 

4.7 X 10-2 

5.1 X 10-5 

3.1’’ 

0 

0 


“Comparison to the data of IStevenson et alj I2Q10I) 
’’Comparison to the data of ILanotte et all I2Q14I) 



Mixing Ratio 


Figure 6. Modeled compositions of a helium atmosphere on GJ 436 b. The 
mixing ratios of important species are shown as a function of pressure. Each 
color corresponds to one species, and each panel con'esponds to a model 
shown in Figure[5]and summarized in TableQ] The panel label indicates the 
color of the synthesized spectrum in Figure]^ In the third panel (“b”), H., 
and CO lines overlap because their mixing ratios are very close. 


IStevenson et alJ (1201 Ol) corresponds to a brightness tem¬ 
perature greater than 1100 K. With our model atmosphere, 
this thermal emission emerges from the deep atmosphere 
(Figure |5]l, which appears to be only possible for a clear 
helium-dominated atmosphere highly depleted in hydrogen 
and having close-to-solar abundances of carbon and oxygen. 
For such an atmosphere, the mixing ratio of methane is 
extremely low, and there is essentially little opacity at 3.6 
pm, thereby creating an infrared window into the deep atmo¬ 
sphere. An internal heat flux ~ 10 times Earth’s geothermal 
flux may thus drive a deep adiabatic temperature profile and 
lead to a high emission flux in this window region. In our 
best-fit model, the emission of a clear helium-dominated 
atmosphere having a very low amount of methane can 
show high brightness temperatures between 3.6 and 4.0 pm, 
ranging from 1100 to 1400 K (Figure |5]l. Comparing with 
the temperature-pressure profile, this corresponds to emitting 
from as deep as 1-100 bars at these wavelengths. Probing 
such deep atmospheres is probably one of the unique features 
of the proposed helium-dominated atmosphere. For this 
reason, a small energy flux from the interior that drives a 
deep adiabatic temperature profile could lead to appreciable 
changes at specific wavelengths of the emerging infrared 
spectrum. 

We note that the grey-atmosphere approximation does not 



Figure 7. Constraints on the composition of the atmosphere on GJ 436 b. 
The color contours show the goodness of fit by comparing atmospheric mod¬ 
els to the dayside emission JStevenson et al.ll2QlQi:ILanotte et alJl2014D . The 
goodness of fit is defined as the jn, where « = 6 is the number of obser¬ 
vations. The axes define the mixing ratios of hydrogen, carbon, and oxygen 
with respect to the total number of atoms, and the mixing ratio of helium is 
Xhc = 1 -Xy{-Xc -Xq. The axis for the carbon to oxygen ratio is shown 
in terms of log(l —Xq/Xq) to emphasize the sensitivity when the carbon to 
oxygen r atio approaches 1 . The u pper two panels compare with the reported 
values of IStevenson et an d the b ottom two panels compare with 

the reported values of lLmiotte et all 1201411 . The parameters not shown are 
marginalized, including the eddy diffusion coefficient ranging from 10 ^ to 
l(f cm^ s’, which are reasonable values for dee p atmospheres according t o 
the free-convection and mixing-length theories IGierasch & Conr^ [198^. 
Both datasets allow H 2 -poor scenarios, but the data of lLanotteTT^ T20T4f 
also allow Hj-rich metal-rich scenarios. 


capture this infrared window, which makes our calculation of 
the temperature profile no longer consistent with the atmo¬ 
spheric composition. As a plausibility check, we calculate the 
total infrared emission flux from the planet, with and without 
the spikes at 3.6-4.0 fim. We And the case with a deep adiabat 
emits 6.6% more energy through this window than the case 
without a deep adiabat. This additional emission energy, if 
solely attributed to the internal heat flux, would correspond to 
an intrinsic temperature of 370 K, much higher than the 60 K 
assumed in the grey-atmosphere approximation. A non-grey 
radiative-convective model is required to accurately calculate 
the temperature profile and the internal heat flux indicated by 
the emission feature. 
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Figure 8 . Com parison of model transmission spectra of GJ 436 b 
with observations jBallard et al.||2010|; IPont et al]l200^ lAlonso et al.ll2008t 

ICaceres et alJI200^ : IKnutson et alJ2011112014^ . The upper panel shows the 

transmission spectra from 0.3 to 20 /im and the lower panel provides a zoom- 
in view for the wave length of 1-2 fim featuring recent HST observations 
IKnutson et alj|2Dl4al) . The model spectra are based on self-consistent at¬ 
mosphere models with parameters tabulated in Table[T] i. e., the red lines cor¬ 
responding to the best-fit model to lStevenson et al.l 201^ and the blues lines 
corresponding to the best-fit model to lLmtotte^T^ 120141) . The dashed lines 
show the transmission spectra of a clear atmosphere, dominated by molecular 
absorption features. The solid lines show the transmission spectra of an at¬ 
mosphere that has an aerosol layer at pressures between 1 and 100 mbar, with 
an aerosol opacity of 0.001 cm^ g“' at 3 /rm. The aerosol optical property 
is computed for a material with a refractive index of 1.001, and a lognormal 
size distribution with a mean dia meter of 0.1 fim and a size dispersion factor 
of 2, using the method detailed in IHu et aP <20131) . The vertical optical depth 
produced by such an aerosol layer is maximally 0.1 for wavelengths longer 
than 3 fim, yielding no impact on the interpretation of the thermal emission 
spectrum. A helium-dominated atmosphere of GJ 436 b must have an aerosol 
layer to be consistent with the featureless transmission spectrum at 1-2 fim. 


One might ask whether probing such deep atmosphere is 
ever possible. We show that this is possible for a helium- 
dominated and hydrogen-poor atmosphere in terms of molec¬ 
ular absorption, because the common opacity sources at 
this wavelength range, including CH 4 and H 2 -H 2 collision- 
induced absorption, disappear in the helium-dominated atmo¬ 
sphere. H 2 -He collision-induced absorption is insignificant. 
The lack of absorbers thus provides an infrared window. One 
might also wonder, given that the model has a lOx solar metal- 


licity, whether clouds c ould form above ^100 bar and block 
this emission window (iBurrows & SharDlll999h . Comparing 
the modeled temperature profiles to the condensation temper¬ 
atures of the materials commonly suggested to form clouds 
in the atmospheres of irradiated exoplanets, including H 2 O, 
Mg 2 Si 04 , MgSi 03 , and Fe, we find that GJ 436 b is too hot to 
have water clou ds, and too cold to hav e silicate or iron clouds 
above 100 bars (ISudarskv et al.ll200^ . A “cloud deck”, if ex¬ 
isting, is then likely to be deeper than what we are concerned 
with. A cautionary note is that strong tidal he ating may raise 
the te mperature of the deep convective layer (lAgilndez et alJ 
l2014h . and raise the altitude of this potential cloud deck, caus¬ 
ing the clouds to affect the infrared spectra. 

4.3. Formation of a Helium Atmosphere 
4.3.1. Current Escape 

GJ 436 b is undergoing transonic hydrodynamic escape. 
The total EUV flux received by GJ 436 b is 7.3 x 10'^ W esti¬ 
mated from Equations (O, and this estimate is fully consistent 
with X-ray and Ly-a flux measurements of GJ 436 (Appendix 
0 . Assuming 7 = 5/3 for an atomic flow, Ccdc = 5x10 m^, 

Kn,„ ~ 1 (valid when the heat is primarily absorbed over a 
broad range of radius below r*), and r* ~ ro, we estimate the 
critical heating rate to be Qc = 6.5 x 10'^ W for a helium- 
dominated atmosphere, and 2.3 x 10^"^ W for a hydrogen- 
dominated atmosphere from Equation (| 6 l). We find that as 
long as 77 > 0.09, the heating rate would be greater than the 
critical rate and the current hydrodynamic escape would be 
transonic. Since the earlier EUV fluxes have been greater than 
the current flux, hydrodynamic escape on GJ 436 b would 
have always been transonic. 

Our calculation of the escape rate is consistent with the re¬ 
cent HST obs ervation of the GJ 436 system in the Lyman- 
a absorption (iKulow et al.ll2014li . The escape rate of H de¬ 
rived from the Ly-a transit light curve is 4 x 10® - 2 X 10’ 
g s“*, assuming an ionization fraction of 0 . 1 . In compari¬ 
son, the standard energy-limited escape rate of a hydrogen- 
dominated atmosphere calculated from Equation ([l]) would 
be 10 ^ ^ 10 ^® g s“' for an efficiency ranging from 0.01 to 1 . 
Even for a very low efficiency, the standard formula overesti¬ 
mates the escape rate. In our model, the current escape rate is 
7.3 X 10^ ~ 7.3 X 10® g s-i for H and 6.8 x 10^ - 6.9 x 10^ 
g s“' for He, for Xu ranging between 10“^ and 10“’, con¬ 
sistent with the observation. These calculations assume a 
heating efficiency of 10 %, but varying the efficiency from 
a few % to 50% does not change the results significantly, 
because the escape rate is limited by Qc for the transonic 
flow. Here, the apparent agreement is only a proof of con¬ 
cept rather than observational confirmation, because the re¬ 
ported absorption depth is measured after the optical transit, 
and ste llar variabilities could signiflcantly affect the pl anetary 
signal (iLoyd & FrancellT014tlLlama & Shkolnikl2015h . Also, 
the reported high yelocity of the escaping hydrogen corona 
(60-120 km s“^) does not indicate the outflow to be super¬ 
sonic by itself. The absorption at high Doppler shifts may 
yer y well be due to absorption in the wings of the line pro - 

flle dBen-Jaffel & Sona Hosseini 2010l:lKoskinen et al.ll2010h . 

or produced by charge exchange with stellar wind protons 
(iHolmstrom et alJl2008HKislyakoya et aTll2014l) . 

4.3.2. Formation Conditions 

We And that the fractionation factor between helium and 
hydrogen {x 2 ) is ~ 0.2 at present for a hydrogen abundance 
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Age [Gyr] 



Figure 9. Fractionation between helium and hydrogen via hydrodynamic 
escape of GJ 436 b, as a function of the present-day Xh and the present- 
day mass of the atmosphere. Panel (A) shows the history of the hydrogen 
abundance of the atmosphere and the fractionation ratio between helium and 
hydrogen by escape (jr 2 ), for an atmosphere presently having a mass of 10“^ 
(black lines) or 3 x 10“^ (blue lines) planetary mass and various present-day 
Xh distinguished by line types. The total mass loss is approximately 10“^ 
planetary mass in these models. Panel (B) shows the calculated initial Xh by 
color contours, assuming an escape efficiency of 10% and an age of 10 billion 
years. We have assumed the atmosphere to be isothermal with a temperature 
of 2000 K. Even a temperature of 1000 K yields only a very minor alteration, 
indicating the isothermal assumption is adequate. We caution that the inte¬ 
rior and the atmosphere may cool down and contract significantly during the 
first billion years, and our calculations may not capture this early phase of 
evolution. When the mass of the atmosphere is 10“^ planetary mass or less, 
hydrodynamic escape can reduce the hydrogen abundance of the atmosphere 
by orders of magnitude. 


ranging from 0.1 to 10“"^, and has remained below ~ 0.2 for 
the entire evolution history (Figure|9]l. The fractionation fac¬ 
tor is significantly smaller than unity, meaning that hydrody¬ 
namic escape disproportionally remove hydrogen from the at¬ 
mosphere. 

The cumulative effect of the selective loss of hydrogen via 
hydrodynamic escape depends on the mass of the atmosphere. 
If the initial mass of the atmosphere is 10“^ planetary mass 
or less, hydrodynamic escape would decrease the hydrogen 
abundance of the atmosphere by more than one order of mag¬ 
nitude on a time scale of a few billion years (Figure 0. The 


change of the hydrogen abundance is greater for a smaller ini¬ 
tial mass. Therefore, the timescale required to significantly 
reduce the hydrogen abundance from around the cosmic value 
to 10“^ ~ 10“"^ by hydrodynamic escape is sensitive to the 
mass of the atmosphere. For example, hydrodynamic escape 
for 10 billion years can fractionate an atmosphere that is 10“^ 
planetary mass or less (Figure |9]l. Similarly, hydrodynamic 
escape for 5 billion years can fractionate an atmosphere that 
is 3 X 10“"^ planetary mass. 

The mass fractionation can be more efficient if the ho¬ 
mopause is well below the EUV absorption altitude and has 
a temperature close to the planet’s equilibrium temperature. 
For a temperature of 600 K, additional calculations indicate a 
fractionation factor between helium and hydrogen of ^ 0.1. In 
that case, hydrodynamic escape for 10 billion years can well 
fractionate an atmosphere that is 3 x 10“^ planetary mass. 

The amount of hydrogen and helium accreted by the 
planet when it formed depends the disk dissipation time, 
the orbital location, and the mass of the planet’s core 

(lRafikov| |2006HHansen & Murravll2012Hlkoma & Horill2012l 

llnamdar & Schlichting|l201^ . If the planet formed at large 

orbital separation and migrated to close-in orbits, iRafikovI 
(l2006l) shows that a 2Q-M^ core would only accrete 0.1 
Mq of hydrogen and helium at 10 AU if the core has a high 
luminosity and the accretion is fast. This makes the ini¬ 
tial mass of the hydrogen-helium envelope 0.5% the plan¬ 
etary mass. Engulfing planetesimals or radiogenic heating 
would increase the core luminosity and reduce the amount 
of hydrogen and helium accreted. After accretion, the ini¬ 
tial atmosph ere may also be re d uced in mass by ph otoevap- 
oration (e.g. iLopez et al.ll2012l lEortnev et akl 120131) . If the 
planet formed in situ dHansen & Murravll2012h . the amount 
of hydrogen and helium accreted would be on the orders of 
10 “^ ^ 10“^ the core mass, and this initial envelope is typ¬ 
ically reduc ed in mass by giant impacts b y 1 ~ 2 orders of 
magnitudes (llnamdar & Schlichting|l20T4h . Both formation 
scenarios could provide conditions to form the helium atmo¬ 
sphere. 

It is also important to check how long the proposed helium- 
dominated atmosphere would be stable against further tran¬ 
sonic hydrodynamic escape. Apparently the total mass loss 
for 10 billion years is on the order of 10“^ planetary mass. 
Therefore, it is conceivable that a helium-dominated envelope 
of 10“^ planetary mass would be stable for billions of years. 
Eigure |9] also shows evolution scenarios in which the planet 
has already had a helium-dominated atmosphere ~ 5 billion 
years ago, and the helium-dominated atmosphere persists to 
the present day. Therefore, although hydrodynamic escape 
might eventually deplete helium from the planet, the evolu¬ 
tionary phase when helium dominates the outer envelope can 
be fairly long. 

We therefore propose that exoplanet GJ 436 b may have a 
helium-dominated atmosphere that evolved from a primordial 
hydrogen and helium envelope. The proximity of the planet 
to its parent star provides the conditions to maintain transonic 
hydrodynamic escape throughout the evolution history. The 
mass and the size of the planet determine whether the escape 
rate has been close to the diffusion-limited escape rate of hy¬ 
drogen. As a result, the planet has experienced dispropor- 
tional loss of its primordial hydrogen. Some of the primordial 
helium has also been lost during this evolution, but some may 
have remained. 
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5. DISCUSSION: DETECTING EVOLVED EXOPLANET 
ATMOSPHERES VIA TRANSMISSION 

The planetary evolution scenario for GJ 436 b can be a gen¬ 
eral process happening on many Neptune and sub-Neptune- 
sized exoplanets. Due to depletion of hydrogen, the abun¬ 
dances of major carbon- and o xygen-bearing speci es would 
change by orders of magnitude (iHu & Seageiil201^ . For ex¬ 
ample, the most important trace gases in an H 2 -dominated at¬ 
mosphere are CH 4 and H 2 O at equilibrium temperatures lower 
than ^ 1000 K (defined for zero albedo and full heat redis¬ 
tribution). When the planet loses hydrogen, the dominant 
carbon- and oxygen-bearing species change to CO, high-order 
hydrocarbons, CO 2 , and O 2 . 

With the recognition that primordial helium may persist as 
a planet loses hydrogen, we propose that the highly-evolved, 
helium-dominated atmospheres can be detected by measuring 
the transmission spectra of Neptune- and sub-Neptune-sized 
exoplanets. The helium atmosphere can be distinguished 
from a hydrogen atmosphere by its molecular compositions, 
i.e., the lack of CH 4 and the dominance of CO, CO 2 and 
(see Figure EJ. These gases have spectral features in 
the visible and near-infrared wavelengths. The helium at¬ 
mosphere, with a mean molecular mass of ^ 4, has a large 
scale height for the s pectral features to manife st strongly in 
transmission spectra (iMiller-Ricci et al.l l2009h . In fact, in 
comparison with recently acquired transmission spectra of 
Neptune- and sub-Neptune-sized exoplanets re.g. iBean et al. 
201T: iBerta et al.l I20i 2|: |&eidberg et al.l l2014t ll^aine et af 
2014 jKnutson et al.ll2014alll^ . observing helium atmospheres 
is within the reach of current observatories, amid complica¬ 
tion of clouds and aerosols masking the molecu lar features 
(e.g. lBenneke & Seageiil201^lMorlev et alJl2013h . 

In addition to transmission spectroscopy, could helium- 
dominated atmospheres be directly detected? We suggest 
the following three methods, taking advantage of the unique 
physical properties of helium. First, helium has a very low 
heat capacity compared to hydrogen, because a helium atom 
has only three degrees of freedom and a hydrogen molecule 
can have six. For example, at 1000 K helium has a specific 
heat capacity of 5.2 J g“' K“’, compared to the hydrogen’s of 
15.1 J mol”' K“' (http;//webbook.nist.gov/chemistry/). This 
would lead to a helium-dominated atmosphere having the 
greatest adiabatic temperature-pressure gradient among com¬ 
mon planetary atmospheres, given that the temperature gradi¬ 
ent of a convective atmosphere is inversely proportional to the 
specific heat capacity. High-resolution infrared spectroscopy 
in the future could potentially measure the temperature gra¬ 
dient by deep spectral features, and thereby verify the dom¬ 
inance of helium. Second, due to its low heat capacity, a 
helium-dominated atmosphere is more likely to have an un¬ 
shifted hot spot, i.e., a hottest zone at the substellar point, 
compared with a helium-dominated one. This will result in a 
thermal emission phase curve peaked at occultation. Third, 
helium is a very poor Rayleigh scatterer, with a Rayleigh 
scattering cross section 20 times smaller than that of H 2 
(iTarafdar & Vardvall969h . As a result, reflection of the stellar 
light by a helium-dominated atmosphere should be dominated 
by multiple scattering of clouds or aerosols in the atmosphere. 
This condition would lack the str ong, broad linear polariza- 
tion peak at a 90° phase angle (iMadhusudhan & Burrows! 

* The transmission spectra of helium atmospheres can still have prominent 
H 2 O features because H -,0 absorption is so strong that a small mixing ratio 


IM3)- 

6 . CONCLUSION 

This paper proposes a new concept that some Neptune¬ 
sized and sub-Neptune-sized exoplanets may have fossil 
helium-dominated atmospheres. This concept is important 
because sub-Neptune-sized planets are ubiquitous in our in¬ 
terstellar neighborhood, and because helium atmospheres 
have unique molecular compositions and can be detected via 
transmission spectroscopy due to an extended scale hight. 

The proposed helium atmospheres can be formed by atmo¬ 
spheric hydrodynamic escape driven by stellar irradiation. For 
hydrodynamic escape to produce a large impact on the overall 
atmospheric composition, the mass of the initial atmosphere 
must be 10 “^ planetary mass or less, and the escape rate must 
be close to the diffusion-limited escape rate. How restrictive 
the first condition is uncertain and depends on the formation 
processes of Neptune- and sub-Neptune-sized planets. We 
show that the second condition is likely to be met by many 
irradiated Neptune- and sub-Neptune-sized exoplanets. 

Applying the concept to the Neptune-sized planet GJ 436 b, 
we find that a helium atmosphere can be formed on a time 
scale of its age, if it has a small atmosphere envelope. He¬ 
lium atmosphere models appear to provide better consistency 
with the planet’s spectra than hydrogen atmosphere models, 
naturally explaining the lack of CH 4 absorption and the dom¬ 
inance of CO as the carbon-bearing species indicated by the 
thermal emission. 

The evolutionary scenarios presented in this paper are or¬ 
ders of magnitude in nature. Refinements are warranted for 
future work. In particular, we neglect the thermal evolution 
and the core contraction in the model, and we assume an 
isothermal atmospheric temperature profile for ease of cal¬ 
culation. These facto rs should be treate d in full evolution¬ 
ary simulations (e.g. iFopez et'd] |20T3)- We also neglect 
atmosphere-surface exchange. If the surface emits hydrogen 
at high fluxes, helium atmospheres cannot form. Therefore 
the evolution scenarios should also be coupled with geophys¬ 
ical calculations to fully map the potential outcomes of the 
evolution of short-period, low-mass exoplanets. Despite these 
caveats, we suggest that observations of Neptune- and sub- 
Neptune-sized exoplanets in the near future would test the at¬ 
mospheric evolution theories presented here, and provide im¬ 
portant clues on the origins of these alien worlds. 
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can lead to large features. 
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APPENDIX 


A BROAD-BAND EMISSION SPECTRUM OF GJ 436 B OBTAINED BY SPITZER 

A broad-band emission spectrum of GJ 436 b in 3-30 fj,m was obtained by measuring the eclipse depths as the planet occults 
its parent star by Spitzer (IStevenson et al.ll2010h . R ecently, the same ob servations were re-analyzed by an independent group, 
and new values of the eclipse depths were reported (iLanotte et al.ll201^ . The two datasets are consistent in the detection of an 
emission flux at 3.6 /rm, and the non-detection at 4.5 pm. However, the specific values and uncertainties for the emission fluxes 
are different between the two datasets. The more recent dataset has a significantly smaller eclipse depth at 3.6 pm, and slightly 
smaller eclipse depths at 5.8 pm and 8.0 pm. 

One should exercise caution when interpreting the reported values, as the multiple-wavelength observations of Neptune-sized 
exoplanets have only started and our understandings of instrumental systematic errors and our data analysis technique are still 
improving. The fidelity of the data may be best established by repeating the measurements and verifying that the results are 
consistent within the quoted errors. For GJ 436 b, this has only been done for the 4.5-pm and 8.0-/im bands. The 8 -/im emission 
flux has been re-observed multiple times an d the results of these repeated observations are consistent with each other at Ict level 
dStevenson et al.ll2010tlKnutson et aTll201 ll) . 

Unfortunately, the measurements of the eclipse depth at 3.6 pm have not converged. The more recent analysis claimed that 
the 3.6-pm light curve of GJ 436 of IStevenson et al.l (120101) exhibits significant variability out-of-eclipse, which may be due to 
incomplete removal of systematic errors. The more recent analysis appears to be able to remove this variabili ty, but derive a muc h 
smaller eclipse depth. However, this analysis used ordinary aperture photometry assuming small aperture sjLanotte,et^|20W ). 
It is therefore unclear how well the intrapixel sensitivity variations have been corrected in their analysis (IStevenson et al.ll2012l) . 
Repeating the measurements of the eclipse depth of GJ 436 b at 3.6 pm would be necessary to pinpoint its emission flux. 


MODEL-OB SERVATION DISCREPANCY AT 16-MICRON 

A remaining challenge for GJ 436 b is the high thermal emission flux of the planet detected at 16 pm, implying little CO 2 
absorption in the atmosphere. Our best-fit models differ from the observation at this wavelength by 2.5 a, and this is the dominant 
source of for the best-fit models (Figure|5]l. We have made several attempts to improve the fit. 

First, making the carbon to oxygen ratio unity or very close to unity does not solve the problem. Although the amount of CO 2 
can be decreased, the amount of CH 4 will be increased for a carbon to oxygen ratio closer to unity. As the 3.6-pm measurement 
for CH 4 has a much smaller error bar than the 16-pm measurement for CO 2 , our fitting still favors a solution that primarily 
satisfies the 3.6-pm constraint. Here we note that further increasing the internal heat flux may su ppress CH,, for large carbon to 
oxygen ratios, and such a large heat flux might be physically plausible if tidal heating is efficient dAgundez et al.ll2M^ . 

Second, we And that modifying the temperature profile will improve the apparent goodness of fit. The results shown in this 
paper are based on grey-atmosphere temperature profiles with the opacities consistent with the composition. The calculation has 
also assumed a dry adiabatic lapse rate of helium for the irradiation-induced convective layer that the thermal emission spectrum 
is probing. The dry adiabatic lapse rate of helium is quite large, which amplifies the size of spectral features including the 
CO 2 absorption feature. As a test we have altered the temperature-pressure profile gradient for pressures smaller t han 0.4 bar, 
mimi cking an adiabatic lapse rate softened by condensable species, or an inefficient day-night heat recirculation (iLewis et al.l 
I 2 OIOI) . We found that the fit to the \6-pm emission flux can be improved to within 2 a without altering the composition. 

EVOLUTION HISTORY OF X-RAY AND EUV FLUXES 

The hydrodynamic escape is driven by of X-ray and EUV radiation from the star. Although the stellar X-ray flux can be 
measured, the EUV flux is not directly observable due to strong interstellar absorption. Eor GJ 436, its EUV flux may be 
estimated either from the X-ray flux or from the Ly-a flux. __ 

The X-ray flux of GJ 436 has been measured to be log(Lx-ray) = 20.16 (iPoppenhaeger et £0120101) or log(Lx-ray) = 18.96 
in the unit of J s“'. Using a fitting relationship between the X-ray flux and the EUV flux 
we estimate an EUV flux of log(LEUv) = 20.13, but this value could vary by 3 orders of magnitude 


Sanz-Forcada et al. 

poiij) 

Sanz-Forcada et al. 

|2onh 


due to uncertainties in the fitting relationship. 

Alternatively, on e could estimate the EUV flux from the L y-a flux. Reconstruction of the Ly-a flux gives an estimate of 
log(LLy-a) = 20.65 (lEhrenrei ch et al.l 201 It iFrance et al.ll2013l) . Using fitting relationships between the fluxes of multiple EUV 
bands and the Ly-a flux (iLinskv et al. 120141) . we estimate log(LEUv) = 20.54, which is consistent with the estimate from the X-ray 
flux. 

Since we are concerned about the evolution of the atmosphere on GJ 436 b, an age-dependent evolutionary history of the X-ray 
and the EUV fluxes is necessary. We adopt the following empirical estimate s for stellar EUV fluxes ba sed on X-ray observations 
and coronal modeling of 80 stars with spectral types spanning from M to E (ISanz-Eorcada et aT]l201 ll) 


log(LEUv) = 22.12-1.241og(r), 
log(Lx-ray) = 21.28-1.44 log(r), 


(Cl) 

(C2) 


where the luminosity has a unit of J s“‘ and the age (r) has a unit of billion year. As a fact check, this adopted evolution history 
gives current fluxes of log(Lx-ray) = 19.84 and log(LEUv) = 20.88, and these values are w ell within the ranges permi tted by the 
above-mentioned observations. The age dependency is also consistent with a recent study (IShkolnik & Barmml2014l) . 
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